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Abstract. We present the results obtained on the magnetisation relaxation 
properties of an XX quantum chain in a transverse magnetic field. We first consider 
an initial thermal kink-like state where half of the chain is initially thcrmalizcd at a 
very high temperature Tb while the remaining half, called the system, is put at a lower 
temperature T s . From this initial state, we derive analytically the Green function 
associated to the dynamical behaviour of the transverse magnetisation. Depending 
on the strength of the magnetic field and on the temperature of the system, different 
regimes are obtained for the magnetic relaxation. In particular, with an initial droplet- 
like state, that is a cold subsystem of finite size in contact at both ends with an 
infinite temperature environnement, we derive analytically the behaviour of the time- 
dependent system magnetisation. 



1. Introduction 

Non-equilibrium properties of quantum chains generated by initial inhomogeneities, 
or by the sudden variation of an external field have been the purpose of several 
extended studies starting from the sixties and continuing up to nowadays. In particular, 
in the context of free-fermionic models, that is Ising-like chains [H [2], since the 
dynamics can be solved exactly [31 S], a lot of exact results have been obtained. 
One may mention magnetisation relaxation profiles which show a scaling behaviour 
13 El CH El El HOI EH EE21 E321 E], two-point correlation functions characterising ageing 
properties [TQl E] , magnetisation and energy current densities [HI [151 EES El EB] > an d 
one may notice that a quantized relaxation behaviour was also observed [T9 | [20 1 121] ■ Since 
the very first studies of Niemejer and Tjion and others [221 ESI EH [25], the relaxation 
obtained shows typically a very slow decay. Indeed, they have shown that the relaxation 
of the magnetisation in a spin-free fermionic chain after a sudden change of the external 
magnetic field is algebraic instead of exponential [26] and that there is a lack of ergodicity 
in the sense that the magnetisation does not relaxe toward its equilibrium value [24J. 
Consequently, in such systems there is no finite relaxation time. Moreover, it was shown 
that for certain inhomogeneously magnetised initial states, at special wave vectors, a 
slowing down of the magnetisation relaxation can occur for gapped free fermionic models 
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[9]. The effect of quenched disorder was also considered in the Ising quantum chain [27J. 
Once more, for an initially completely ordered state, either in the field direction or 
perpendicular to it, an algebraic decay of the magnetisation toward a stationary state 
was obtained. More recently, in the critical Ising quantum chain we have studied the 
relaxation of an initially magnetised finite droplet at a given temperature in contact 
at both ends with infinite Ising chains. We have shown that the droplet magnetisation 
decays in time to the equilibrium magnetisation with a l/tf(x/t) behaviour [20] where 
f(u) is a scaling function depending on the temperature. For a vanishing droplet initial 
temperature, the scaling function is simply given by the characteristic function Tl(u). 

In the present study, we focus our attention to the same situation in the XX 
quantum chain in a transverse magnetic field. Whereas in the Ising case, the transverse 
magnetisation is not conserved in time, in the XX quantum chain it is a conserved 
dynamical quantity and one may expect that this conservation law will show up in 
a way or an other in the relaxation behaviour. To be more specific, let us precise 
the initial set-up. The system is initially prepared in a factorized canonical state 
p(0) = pb(0)p s (0) where the density matrices Pb )S (0) oc e~^ b ' aHh - s are canonical states 
of parts of the (infinite) chain each thermalized at a given temperature j3^ 8 . At time 
t = 0, the (local) interactions between the different parts are switched on and the 
dynamics is generated by the total Hamiltonian H, = Tib + 7i s + Tii. The paper is 
organised as follows: In the next section we give some details of the by now classical 
diagonalisation of the XX-model and its dynamics. Section 3 gives more details on 
the initial set-up and on the time- dependent expectation values of physical quantities, 
putting a special focus on the transverse magnetisation. In particular, it is shown that 
the time-dependent transverse magnetisation may be expressed as a discret convolution 
product of a Green function with the initial magnetisation profile. In section 4, from 
an initial kink-like state we derive analytically the magnetisation Green function and 
show that, at vanishing system temperature, one has to distinguish between the low 
and high field regimes. Consequently, in section 5, we show that the relaxation of the 
magnetisation of an initial finite droplet of size L s , toward its equilibrium vanishing 
value, changes after an initial linear decay from an algebraic decay, L s /t, at large fields 
to a faster size-independent decay at lower fields. However, this behaviour does not 
survive for finite system temperatures, but nevertheless the magnetisation relaxation 
at low enough temperatures is strongly affected. We finally summarize our results and 
draw out some conclusions. 
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2. Diagonalisation and Dynamics 



2.1. Diagonalisation 



The free boundary conditions XY quantum chain in a transverse field is defined through 
the Hamiltonian 

L-l i L 

H = ~ 2 E + VM+i] - 2 E < w 

n=l n=l 

where the as are the Pauli matrices at site n, /i is the transverse field applied in the 
z-direction and where we parametrize the coupling constants using the anisotropy factor 



k G [—1, 1] as J x 



k)/2 and J y 



[1 — k)/2. In the special case k 



the ^-direction coupling constants vanish and one recovers the Ising quantum chain 
case. At K — 0, J x = J y and the Hamiltonian describes the XX-model for which 
the total transverse magnetisation M z = ^ n <5 n is a conserved quantity since we have 
[H, M z ] = 0. Consequently, the exchange interaction energy and the Zeeman energy are 
separatly conserved quantities. 

In terms of the lattice Clifford's operators r* , 

n-1 

T n = H(-<i)< 

n-1 

T n= ~ Ui-^X ( 2 ) 



with rtj 



T l n satisfying the anticommutation rules 
the Hamiltonian (TjQ) takes the quadratic form [28j tH E] 



n 



4 



(3) 
(4) 



where T is the 2L-component vector 



with r 1 




ri 



and i = 1,2. 



(5) 



V n ) 

The 2L x 2L hermitian matrix T is given by 

C 

Ct 



(6) 



with 



/ h 

J, 



\ 



■-hi 

h J„ 



h 

J, 



\ 



h ) 



(7) 
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Parametrizing the eigenvectors V q , associated to the eigenvalue e q , of T as 

we obtain from the eigenvalue equation TV q = e q V q the following equations: 

- iCtpq = e q (f> q , C ] (p q = -ie q i) q . (9) 

We can notice here that if we change simultaneously e q — > — e q and tp q — > — ^ g these 
equations are unchanged. So, to each eigenvalue e g > corresponds an eigenvalue 
e ? / = — e q with the associated eigenvector V q i((j) q i ,ip q >) = V q ((j) q , —ip q )- 

Introducing the fermionic creation and annihilation operators n^, rj q 



^(0 g (n)r^ + # g (n)r 2 n ) (10) 



satisfying the canonical Fermi-Dirac anticommutation rules {m,T) q >} = 5 9i9 /, the 
Hamiltonian ([1]) takes the free-fermionic structure 

L 

« = ZNMta - 1/2) (ii) 

9=1 

where e q are the positive eigenvalue of T, given by 



'K 2 sin 2 (g) + (/i + cos(g)) 2 . (12) 

In the XX-quantum chain, with k — 0, one has for the (positive) spectrum 

e q = \h + cos q\ = \h — cos k\ (13) 

with k = 7T — q £ [0, 7r] and where the gap vanishes for h < 1 at kp = arccos(/z). 
In an equivalent picture, the (positive) spectrum associated to k £ [0, fop] is usually 
interpreted as corresponding to the creation of holes in the Fermi sea, while the ground 
state corresponds to rifc=o^llO) wnere the W are particle fermionic creation operators 
with negative energies = — e^. 

2.2. Dynamics 

The Heisenberg equations of motion for the lattice Clifford's operators are easily solved, 
see the appendix 1 for more details, and lead to [2] 

rut) = e Mt T{e- int = ^nnmi (u) 

k.u 

where (C\D) = |{C',£)} is a pseudo-scalar product with {.,.} the anticommutator. 
The time-dependent contractions are explicitely given by 

(rK(t)>= Y^Uk)Un)^{e q t) 
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(rl\r 2 n (t)) = J2Uk)Mn)™s(e g t) 

q 

(Tl\rl(t)) = (T 2 n \rl)_ t = -J2 MWM sin M • ( 15 ) 

q 

In the XX-chain case, in the thermodynamic limit they take the form 

«l r iW> = ^ J -w{^S); = (16) 

(M) =,-^ w { (in 

with (r||r^(t)) = (r^|r^(t)) and where J n (i) is the J-Bessel function of order n. These 
contractions are the basic time-dependent quantities of the XX-chain. 

Formally, since {Tji.|r{) = Sij5 nt k, the set of {T k } forms an orthonormal basis of 
a 2 L- dimensional linear vector space £ with inner product defined by (.|.) = > •}■ 
Hence, every vector X e £ has a unique expansion X = k (T k \X)r % k . A string 
operator of the form O = XiX 2 . . . X n is a direct product vector in £® n and has the 
expansion 

O = X,X 2 . . . X n = ■ ■ • E K ■ ■ • ( r fe \X n )Tt • • • l£ . (18) 

3. Initial conditions and expectation values 

3.1. Product initial state 

The system is initially prepared in a factorized state (a product measure state) 
PHI H21 031 EDI EI] of the form 

p(o)=n^(°) ( i9 ) 

i 

where the density matrix 

p i (0) = i-oq)(-/3 i W i ) (20) 

is a canonical state, at inverse temperature (3j, associated to the Hamiltonian Tij of the 
jth subsystem of the total chain decomposed in the following way: 

n E n * ' E h 'j ■ • ( 21 ) 

3 3 

where the Pauli's matrices of the j th part are written (o~ n )j. The interface interaction 
term is 

Wii+i = -J [K,);K); + > + K , (22) 

where we wrote the Pauli matrices of the j th part as (<r n )j- 
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3.2. Expectation value of the transverse magnetisation 

The expectation value of an observable O at a time t is given by 

(O(t))=Tr{O(t)l[p J (0)}, (23) 



where G(t) = e lHt Oe lHt is the operator associated to the observable at time t in the 

In ' 



Heisenberg picture. In particular if O is a string of Clifford's operators O = Tj 1 . . . T\ n 



one has at time t 

<°(*)> = E • • • EfltreW) • • • < r fel r t(*)>Mrt • • -r-p(o)} , (24) 

where the elements TrjTj^ . . . T l k n p(0)} represent the relevant initial state information 
(correlations) for the time evolution of the quantity O. Since the time-contractions are 
explicitely known one needs only to compute the initial state properties of the string 
operators for a given p(0). 

In this work we focus our attention on the time behaviour of the transverse 
magnetisation which is, as already stated, a conserved quantity. The transverse 
magnetisation is one of the simplest Clifford string operator. Namely, it takes the 
form = — ir^rjj and at time t, in the Heisenberg picture we have 

<(t) = -K(tK(t) . (25) 

Its expectation value, in the thermodynamic limit, at time t is given by 

(0(*) = E*?(*,*)w ^ 

where h,k+j = Tr{—iY 2 k Y\ + -p{Q)} is the initial condition matrix. The time-dependent 
functions F 3 t are given by 

F>(n,k) = Fl\rm<fi^K(t)) - (nK(t))(n +J K(t)) ■ (27) 

and they can be explicitely evaluated in the thermodynamic limit: 

F? j (n, k) = F 2 t \n — k) = (-1)' J n - k (t) J n - k - M {t) , (28) 

and they identically vanish for odd j. Consequently, we have for the time behaviour of 
the transverse magnetisation 

oo 

K)(t)= E (n), (29) 

j=-oo 

where * stands for a discret convolution product and where we introduced the notations 

One needs now to evaluate the initial matrix J. First remark that for Clifford 
operators T belonging to different initial blocks, the trace identically vanishes, reflecting 
the fact that they are traceless operators and that there are no initial correlation between 
them since the initial state is of a product form. So, the only non- vanishing contribution 
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comes whenever the two-rs are in the same sub-system. Consequently, the / matrix has 
a diagonal form such that 

/-. \ 

h-i 

Ii 

V 

where l\ is a L\ x L\ matrix with L\ the size of the Zth subsystem, with matrix elements 

'#4 



(30) 



(ii) k>k > = -^24^)4^')^ 

9=1 

where the /-labeled quantities refer to the Zth subsystem Hamiltonian 7ii 



(31) 



4. System in contact with an infinite temperature "bath" 

4-1- Transverse magnetisation Green function 

We start our analyse by considering a chain where half of it (the left) is thermalized at 
inverse temperature fa, and is called the bath, and where the right half, the system, is 
initially at inverse temperature fa = fa Consequently, the initial matrix I is 

I 6 
I, 



I 



(32) 



Taking the infinite bath temperature limit, the initial matrix I& vanishes so that only 
the system part remains with 

(/.)*,* = - E m k W k ') tanh (x) (33) 

where L s is the system size. The system initial matrix can be decomposed into a surface 
contribution, near the interface boundary, and a homogeneous bulk part: 



with 



s)k,k' 



Bulk\ 



(I 



Bulk\ 



lk,k' 



a« Surf )k,k' 



n=X 



cos (q n (k - k')) tanh {j3e q J 2) 



[l s Surf ) k y 



E l + T cos ^ + kl)) tanh (/5e ""/ 2) • 



In the large size limit, L s 



n=l 
OO, 



(I 



Bulk\ 



)k,k' 



\k+k' 



w de 



71 



cos ((k- k') 9) tanh (J3e h (6)/2) 



and 



fi s 5ur/ : 



— cos ((Jfe + jfc')0) tanh (j3e h (0)/2) 



7T 



(34) 

(35) 
(36) 

(37) 
(38) 



for the surface contribution, with e/ l (6 l ) = h — cos(6). 
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4-2. Zero temperature limit 

In the system low temperature limit, T s = 0, one has to distinguish between the region 
h < 1 and h > 1 since tanh(/3e /l (^)/2) = sgn(e /l (^)). 



i) /i > 1 region For h > 1, the excitations are always positive and the sum in (1331) 
reduces to — ^2 ipg(k)(p s q (k') = 5k t k', so that the initial matrix is just the identity matrix. 
The average magnetisation, in the thermodynamic limit L s j> — > 00 is then given at later 
time by 

00 

m*(n,t) = = F t( n - k)H(k) , (39) 

k=— 00 

where H is the Heaviside function. The Green function G t {n) = F^{n) = J 2 {t) can be 
expressed in the continuum limit as p2] Gtin) ~ \g{n/t) with the scaling function 

= /T~ 2 ' l U l ^ 1 ^ 

7TVl — U 

and zero otherwise. Finally, with the step-like initial condition we have in the continuum 
limit 

{0 x<-t 
^ arcsin(x/t) + \ —t<x<t . 
1 t <x. 

This solves the extrem temperature set-up for h > 1. 

ii) h < 1 region For ft, < 1 one has a ^-independent profile for the bulk contribution 



where # F = arccos(/i), and for the diagonal part, j = 0, 

p° = m 2 (0) ee (I s Bil ' fc ) fc>fe = - arcsin(/i) . (42) 
The surface term leads to a fc-decaying part 



(1 Surf, f 1 , j 2 sin((2fc + j)g F ) 



Wi=W- ^— Vj. (43) 

Neglecting the boundary effect term, I Sur f j and using the fact that if i+1 (n) = 0, one 
has for the transverse magnetisation 

' * m z (0)H I (n) = (G t * m 2 (0)) (n) , (44) 



«>(0 



^m z (0) * 



where we have used the fact that (l Bulk )k,k+2j = P 2j (k) = p 2j H(k). The Green 
function Gt is determined by the spatial derivative of the former equation, since 
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Figure 1. Magnetization profiles with T s = 0, T& = oo for different fields h = 1. 
h = 0.8, h = 0.6 and h = 0.4 from top to bottom on a chain of total size L = 500. The 
analytical results are in full lines. 



[m z (x,t)]' = G t * m z (0)H' = m z (0)G t * 5 = m z (0)G t . In the asymptotic limit (see 
the appendix) one obtains 

G t {x) = -g (-) = 1 , < |-| < 1 , (45) 

V; r\tJ t2arcsin(/ i ) v /l-(x/t) 2 ~ l t l ~ ' v ; 



and zero otherwizeJ§| This disconnected support of the Green function will play an 
important role when considering the relaxation of a finite initially magnetised domain. 
It will lead to the appearance of a finite relaxation time. 

It turns out that the step-like initial state generated by the two-temperature set-up 
develops a plateau of constant magnetisation, as it was already observed in ([161 fl3])- 
due to the vanishing of the Green function in the domain \x/t\ < y/l — h 2 : 

m z (x, t) = $ h {x/t) + m z (0)/2 (46) 

with $h(— u) = —&h(u) where 



< u < VI - h 2 

$h{u) — { m z (0)/2 - arccos(n)/7r Vl - h 2 < u < 1 . 

m z (0)/2 \<u 

We show the numerical and analytical results in figure ([1]) for different fields. 
4-3. Finite temperature effects 

While keeping a very hot bath (formally taken to be at infinite temperature) we now 
consider the effect of a finite initial temperature system state on the relaxation of the 

§ Note that the same analysis is handablc in the Ising case at critical transverse field. The result 
is simpler with a Green function given by Gt{x) = jg Ism s = J_n (^-) where IT(u) = 1 Vit € 
[—1/2, 1/2] and zero otherwize. Details of the derivation are given in appendix 3. 




Figure 2. Numerical rescaled magnetisation profiles for three different inverse 
temperatures (At x/t = from top to bottom /3 S = 1, (3 S = 10, j3 a = 100) at h = 0.8. 
The size of the system is L s = 100 and the time is t = 700. The dashed lines are the 
corresponding analytical curves for the Green functions. 



transverse magnetisation. As in the previous section, we can compute the Green function 
in the asymptotic limit from the step-like initial condition. One finds for h > 1 (see the 
appendix for details) after a lengthy derivation 



G t (x) = 

27rtm 2 (0) 




(47) 

where j3 s is the inverse initial system temperature and where 

m z (0) = I jf d9 tanh (h - cos 0) J (48) 

is the initial bulk magnetisation density. Expression (|47|) is also valid at h < 1. On 
figure (j2J) we have plotted the numerical magnetisation profiles compared with analytical 
result (1471) for different temperatures, since in the long time limit t ^> L s one recovers 
the Green function as m z (x,t) = G t (x) * m ^ Ls 8 {x/t). 

One may notice that it is also possible to deduce the Green function in a very 
simple way starting with the kink-like initial state and assuming that the magnetisation 
at the interface takes for all times its asymptotic value m z (0)/2, since the initial bath 
magnetisation is zero while the initial system magnetisation is m z (0). Indeed, using 
( I44jl . which is in fact also valid at h > 1, one has m z (x,t) = (Gt*m z (0)H)(x). Together 
with the scaling assumption G t (x) = tm l^ g(x/t) one has 

m z (x,t) = J du g(u) . (49) 



Relaxation in the XX quantum chain 



11 



Now, taking into account the fact that at the interface x = 0, the magnetisation is given 
by fj48|) . one has 

m z (0,t)= [° dug{u) = ^-^- = [* —tanh(—(h- cos 9)] . (50) 
J-i 2 J 2tx V 2 / 

Decomposing the integral into two peaces, one from to n/2 and the other one from 

7r/2 to 7r, by a simple change of the integration variable one obtains 

<-° '° tanh(f (h - VT^u 1 )) + tanh (f (h + VT^u 1 )) 

which gives by identification g(u) and so G t (x). 



dug{u)= du ^ " v 2 ^(51) 

7-1 27TV 1 — u 



5. Relaxation of an initial droplet 

If one considers initially a finite system of size L s (the initial droplet) at an inverse 
temperature j3 a in contact at both ends with infinite temperature "baths", from the 
convolution of the asymptotic Green function with P°(k) = m z (0)n one obtains 
the transverse magnetisation at later times. The total droplet magnetisation density 
mZ {t) = £- SfeeL s mZ (k, t) will relaxe toward zero since the baths are at infinite 
temperature (and consequently have a vanishing transverse magnetisation). From the 
asymptotic Green function obtained in the different regimes h > 1 and h < 1 we obtain 
two different relaxation behaviours. 



5.1. Zero temperature droplet 

In the h > 1 regime the droplet magnetisation density behaves as 

l-ff t<n 

m z {t) = { o ', Ls , 2t A A 7ITi\ , (52) 



- arcsm 



with n = Asymptotically the magnetisation relaxes with a power law behaviour 

ma 

m z (t) ~ — h>l. (53) 

One may notice that in the critical Ising chain one has exactly the same relaxation 

behaviour for the droplet mag netization: m z (t) ~ h [2H1I2T] . 

On the contrary, in the low field regime, h < 1, the time behaviour of the droplet 

magnetisation density obtained from the asymptotic Green function ()45l) is given by 

2 / \ 2 t 

m z (t) = — arccos ( vl — h 2 ) — h, t < T\ 

71 \ J TV L s 



m z (t) = — arccos (^/l — h 2 ^j arccos 




n < t < r (54) 




Figure 3. Relaxation of the magnetisation density at fields, from top to bottom, 
h = 1 and h = 0.8 for an initial droplet of size L s = 50. The full lines correspond to 
the analytical results. 



and 

m z (t) =0, t>r 2 . (55) 



with T2 = L s /yl — h 2 = L s j cos(m z (0)7r/2). The vanishing of the magnetisation at 
times t > T2 is a direct consequence of the vanishing of the Green function in the region 
— \/l — h 2 < x/t < \/l — h 2 . However, this exact cancellation is valid only in the limit 
L s — > oo. For a finite droplet, one has finite size corrections to the Green function such 
that at the leading order the Green function is given by G t (x) = g(x/t)/t[l + 0(1/L S )]. 
These corrections lead for the magnetisation density to a size independant correction to 
zero of order \jt. Consequently, the ratio between the magnetization densities at h < 1 
and h > 1 is asymptotically (in time) given by 

"*<'<'> -o(l-). (56) 



m l>i(t) \ L s 

It is in this sense that the magnetisation density within the droplet changes from a slow 
relaxation fl53l) . depending on the size of the system, to a much faster relaxation since 
it is independent of the size. If one considers the total remaining magnetisation of the 
droplet, we pass from a diffusive like behaviour L 2 /t at h > 1 to a balistic-like behaviour 
L s /t at h < 1. These two regimes reflect the fact that the dispersion relation changes 
from a parabolic shape, e(k) ~ h — 1 + k 2 /2 at h > 1 to a linear regime e(k) ~ \k — kp\ 
for h <l. 

On figure ([3]) we show the numerical results compared to the analytical expressions 
obtained for the droplet magnetisation density. The agreement between the exact 
numerical results and the analytical asymptotic expressions is excellent. 
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5.2. Finite temperature droplet 



With the help of the finite temperature Green function ( 1471) . one may compute the 
droplet magnetisation density at the leading order: 



m z (t) 



LA 



( tanh 
7rt 2 V 



0s 



(h 



tanh 



0s 



(h 



(57) 



which shows the L s /t power law decay for any transverse field h. In the h < 1 region, 
the faster relaxation is recovered in the low temperature limit S — > oo. Indeed, for 
sufficiently small temperature the magnetisation density is given by 

m *(t)~^exp(-Ps(l-h)) (58) 

Tit 

and one may identify a typical size-dependent inverse temperature 0(L S ) = In L s /(1 — h) 
and one recovers the fast relaxation for inverse temperatures S ^> 0(L S ). 



6. Summary and discussion 

We have presented here analytical results, confirmed by exact numerical calculations, of 
the relaxation of the transverse magnetisation of the XX quantum chain. In particular, 
starting from an initial product state depending on the strength of the initial magnetic 
field, we have obtained the asymptotic Green function associated to the transverse 
magnetisation. When the system is divided initially into two parts, a hot bath and 
a cold sub-system, we have shown that the Green function takes the scaling form 
Gt(x) = jg(x/t), where the explicit shape of the scaling function g{u) depends on 
the initial temperatures and on the field. In the extreme scenario, with a vanishing sub- 
system temperature and an infinite bath temperature, at low fields, h < 1, the scaling 
function shows the very peculiarity that its support is disconnected. This property leads 
to the emergence of a flat magnetization profile when one starts with a kink-like initial 
profile [TB]. In the case of an initial finite magnetised domain it leads to the fact that its 
relaxation is much faster in the low field regime, h < 1 than it is at higher fields. Namely, 
one has m h< i(t) / m h >i(t) ~ l/L s where m{t) is the domain magnetisation density at 
time t and where L s is the size of the domain. 

Finally, one may notice that for a Green function taking the scaling form Gt{x) = 
jg(x/t), one has identically 

d t G + d x jG = (59) 

so that the impulsional response magnetisation mf = m z (0)G t decays in time and satisfy 
the continuity equation 

d t m z i (x,t) + d x j i (x,t) = (60) 

with the current ji simply given by ji(x,t) = jm z (x,t). In the step-like situation, the 
scaling form jg[x/t) of the Green function leads to a scaling solution m*(x,t) = &(x/t) 
which obviously satisfy the homogeneous scaling equation 

{td t + xd x )m z s {x,t) = (61) 
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which leads for the current j s to the equation d x j s = jd x m z s . So that the step-like 
current is related to the impulsional current through d x j s = j l since d x m z s = mf. 
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Appendix 1: Dynamics of the lattice operators 

The Clifford operators , can be expressed in terms of fermionic operators. 

L 

ri(f) = X>,(n)(#) + 

9=1 
L 

r£(*)= £Wn)(#) -»*,(*)) ( 62 ) 

q=l 

The time evoution of the creation and annihilation operators is obtained by m(t) = 

uUt)4u q (t) with 



U g (t) = e- ie ^ (63) 



q 

and leads to 

Vq (t) = . (64) 

With the help of these relations we can re-expressed the Cliffords operators in terms of 
the initial time operator T. Finally one obtains 

rat) = j^innhn (65) 



with 



9 

(1*1^)*= E^W^Hcos(e^) 
q 

(Tl\r 2 n ) t = (Tl\rl)_ t = -^MWMsH^t) ■ (66) 
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Appendix 2: Zero temperature XX Green function 

The contribution of the bulk to the magnetisation is given by: 
m>,t)=^^F/(n-£;)(I s 

j k 
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Bulk\ 
L s }k,k+j 



(67) 



Using {l s Bulk ) Kk+J = H(k)pi and the relation (I s Bulk ) k , k+J = (l s Bulk ) k+hk one has: 

oo oo 

m z (n,t) = m z (0) ]T J 2 n _ k (t) + 2 ]T (-1)' J n _ k (t) J n - k - 2j (t)p 2j . 

k=0 k=0,j=l 

with p° = m 2 (0). The spatial derivative + 1, t) = m z (n + 1, t) — m z (n, t) is 

(oo ^ 
m z (0)J n (t) +2j2(-iy P 2 iJ n _ 2j (t) 

Using the asymptotic expansion of the Bessel function for v » 1 

v 



(68) 



COS (3 



7Ti/ tan /5 



cos ip 



(69) 



(70) 



with if; = z/(tan/5 — /3) — 7r/4, and taking n — > oo one has 

2 



$'(ra,t) 



7Tt A /l-(f) 2 



X 



COS 2 ('0) 



(0) + 2 ^(-l)V j cos(2j arccos(-)) 



— cos(-0) sin( , 0) 



(71) 



2 ^(-l)V j sin(2j arccos(-)) 
. j>i 

Keeping only the first term in this expression, where cos 2 (^) is replaced by its average 
value 1/2 and taking the explicit expressions for the profiles p 7 , given in ()4T|) . together 
with 

V7 sm(jx) 



one obtains finally 



y^(— 1) J = -x/2 - 7r < a; < 7r , 
g^sin^x) = ^ _ ^ < x < 3?r 



(72) 



(73) 



&(n,t) 



jt < VT^W 



ntVT^} 1 VT^W < n/t < 1. 



for positive n and $'(— n, t) = $'(n, t). The Green function is simply given by 
G(z,t) = &(x,t)/m z (0). 
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Appendix 3: Critical Ising Green function 

In the Ising case the bulk contribution is given by (see [20, 12T]): 

m '( n > *)=EE ^■* +J "(*)(lf H, *)*,*+j, (74) 

with 

(-IJipfc-M-i^) = J 2(n _ k) (2t)J 2{n . k) . 2j (2t) 

— ^2(n-fe)+l(2t) J2( n -fe)-2j-l(2t) . (75) 

At system vanishing temperature T s = the 7 S elements are given by: 

/ t Bulk\ * \ L ! C7a\ 

(Is )fc ' fc ' " *2(lfc'-lfc) + l ' (76) 
Using the relation (I s Bulk ) k , k , = (I 3 Bul %+i, k one has: 

oo 
i=0 fe 

with 

(-l) j Fi(p) = 2 (jj,(2i) J 2p - 2j -i(2t) - J 2p (2t)4_ 2j _ 1 (2t)) . (78) 
The spatial derivative function $'(n + 1, t) = m z [n + 1, t) — m z (n, t) is then given by 

4-\ ^ T/ /ojA ^2n-2j-l(2t) 4 . ^-^ ^2n-2j'-l(2*) /,__,. 

$ n, t) = -J 2n (2t) > — — ^— J 2n (2t) > — — ±— . 79 

7T ' 27 + 1 7T 27 + 1 

Taking this expression together with the asymptotic expansion of the Bessel function 
we finally arrive at 

G(x,t) = — n (-) . (80) 



2t \2t 

Appendix 4: Finite temperature XX Green function 

In the system non zero temperature case, one has 

$'(n,t) = m z (0)J 2 n (t) + J n (t)A n (t) (81) 

with 

oo 

A n (t) = 2j2(-l) j P 2j Jn-2 j (t) (82) 

3=1 

where p 2j H{k) = (I Bulk ) kjk+2 j is defined in (j3"Tj) . Using the relation 

oo 

tanh(x) = 1 + 2 ^(-l)V 2fca; Vx > (83) 



k=l 



(imposing h > 1) one has: 



p 2 ^ = J 2j (0) + 2 ]T(-1) V 2Ah %(/U) (84) 
k=i 
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where l2j{u) are the modified Bessel functions. Then 

oo oo 

A n (t) = 2^(-iy2^(-l) fc e-^ fe %(/5 s fc)J n _ 2j (t) 

3=1 k=l 

The asymptotic behaviour of the Bessel function leads to: 



Jn(t)A n (t) 



wty/i - (n/ty fri 



k -2f3 s hk 



X 



2 ^(-l) j hjiPsk) cos(2j axccos(n/t)). 

3=1 

Since hj(u) = (— 1) J J 2 j(W) one has 



(8 



2 iyi2j(P s k) cos (2j arccos(n/t)) = 2 J 2 j(z/3 s /c) cos (2j arccos( 

= cosh f/3 s AV 1 - {n/t)A - J (ip s k) 
Finally, reinjecting this expression into the previous one, one obtains 

1 



tanh I h 



h 



+ tanh I — 
2 



2 

7 



(8 
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